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ABSTRACT 

Context. IGR 117511-3057 is the second X-ray transient accreting millisecond pulsar discovered by INTEGRAL. It was in outburst 
for about a month from September 13, 2009. 

Aims. We analyze the spectral and timing properties of the object as well as the characteristics of X-ray bursts with the aim to constrain 
the physical processes responsible for the X-ray production in this class of sources. 

Methods. The broad-band spectrum of the persistent emission in the 0.8-300 keV energy band was studied using simultaneous 
INTEGRAL, RXTE and Swift data obtained in September 2009. We also describe in the 2-120 keV energy range the timing properties 
such as the outburst light curve, pulse profile, pulsed fraction, pulsed emission, and time lags, as well as study the properties of X-ray 
bursts discovered by RXTE and INTEGRAL and the recurrence time. 

Results. The broad-band average spectrum is well described by thermal Comptonization with an electron temperature of kT^ ~ 25 
keV, soft seed photons of fcrti, ~ 0.6 keV, and Thomson optical depth Tt ~ 2 in a slab geometry. During the outburst the spectrum 
stays remarkably stable with plasma and soft seed photon temperatures and scattering optical depth being constant within errors. We 
fitted the outburst profile with the exponential model, and using the disk instability model we inferred the outer disk radius to be 
(4.8 - 5.4) X 10'° cm. The INTEGRAL and RXTEAsXa reveal the X-ray pulsation at a period of 4.08 milliseconds up to -120 keV. The 
pulsed fraction is shown to decrease from ~22% at 3 keV to a constant pulsed fraction of ~ 17-1 8% between 7-30 keV, and then to 
decrease again down to ~13% at 60 keV. The nearly sinusoidal pulses show soft lags monotonically increasing with energy to about 
0.2 ms at 10-20 keV similar to those observed in other accreting pulsars. The short burst profiles indicate hydrogen-poor material at 
ignition, which suggests either that the accreted material is hydrogen-deficient, or that the CNO metallicity is up to a factor of 2 times 
solar. However, the variation of burst recurrence time as a function of m (inferred from the X-ray flux) is much smaller than predicted 
by helium-ignition models. 

Key words, pulsars: individual (IGR J1751 1-3057) - stars: neutron - X-ray: binaries - X-ray: bursts 



1. Introduction discovery of a new class of accreting NSs, accreting millisec- 

ond X-ray pulsars (AMXPs). The first direct measurement of 

■ - - > The discovery of radio milUsecond pulsai's in binary systems ^pj^.^p NS during an accretion phase was published by 

m the 1970s lead to the prediction that neutron stars hosted in palanga et al. (2005a) in the source IGR J00291+5934, sti-ength- 

low mass X-ray binary systems (LMXB) were their progen- ening'the hypothesis of the recycling of old r adio pulsars t o 

itors dAlEaiLet al. 1982). As a binary system evolves through mHiisecond periods (see also the reanalysis bv lPatrunollIoToh . 

phases of accretion onto the neutron star (NS), it gains an- Qn the other hand, between the outbursts, long-term monitor- 

gular momentum from the accreted material, which is suffi- jug ^^ows some AMXPs to exhibit spin-down in quiescence 

cient to spin-up the NS to a rotation period equilibrium in ,Hartman et al. 2009; Patruno 2010; Papitto et al. 2010b). The 

the millisecond range. The first confirmation that LMXBs can gpin frequencies of AMXPs lie in the range of 180-600 Hz and 

host rapidly rotating NSs was the discovery of coherent os- orbital periods are between 40 min and 5 hr (see reviews by 

dilations from type-I X-ray bursts (see IStrohmaver & BildstenI |Wiinandsll2006l : lPoutanenl2006l) 

|2006t for a review). Later identification of SAX Jl 808.4-3658 ^MXPs exhibit many characteristics similar to other 

as a 401 Hz pulsar dWijnands & van der Klis| | 1998[) led to the lMXBs. Their broadband specti-a show soft t hermal an d hai'd 

Comptonized components (see reviews by iPoutaneij 120061 ; 
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(iBarret et alJ l2000l) . AMXPs also sh ow quasi periodic oscil- 
lations a nd X-ray bursts (see e.g., 'Chakrabartv et al.l l2003t 
IWatts et al. 2005) . The pulse profiles are typically close to sinu- 
soidal with modulations of 3-15%. The pulse profiles are energy 
dependent a nd demonstrate s oft time lag in the range up to 100- 
250 yus (e.g..lCui et al.lll998l:lFbr d 2000; Poutanen & Gierliiiski 



l2003t iGierliriski & Poutanenll2005c iFalanga & Titarchukii2007i) . 

The magnetic field of the NS as inferred from accre- 
tion m odels is relatively weak, in the range of ~ 10^ - 
10' G (iPsaltis & ChakrabartvllT999t iDi Salvo & Burderill2003h . 
Measurements of the inner disk radius from evolution of the 
pulse profiles and the spin-down of the pulsars between the out- 
bursts allow a more accurate determinati on of the magnetic field 
of 1.5 X 10^ G in SAX J 1808.4-3658 (llbragimov & PoutanenI 
12009^,'Hartm an et al.l2009l) and of 2x 10^ G in IGR J00291-H5934 
(Patruno 201^. 

Pulsations are not detected during quiescence in either X- 
ray or optical bands, nor has milhsecond rad io pulsed emis- 
sion been found (see e.g., Ilacolina et al.l 1201 Ol and references 
therein). Following outbursts, the optical counterpart has been 
observed to diminish in intensity consistent with the X-ray 
flux decay rate implying th e origin of the emis sion to be in 
the accretion disk (see e.g.. iFalanga et al .1 12005 al) . Companion 
stars of AMXPs a re hig hly evolved white or brown dwarfs 
dPelove & Bildstenll2003h . In the AMXPs with an orbital period 
in the hour range and a hydrogen-rich do nor brown dwarf, X- 
ray bursts have been detected (see e.g., iGallowav & CummingI 
12006) . during which burst oscillations are ne arly in phase with 
the coherent oscillation at t he spin period dChakrabartv et al.l 
l2003HStrohmaver et al.ll2003i) . 



1.1. The source IGR J1 751 1-3057 

Among roughly one hundred LMXBs hosting a NS, 
IGR J1751 1-3057 is the twelfth known to host an AMXP 
It is the second source discovered by INTEGRAL during the 
Galactic Bulge monitoring program (Baldovin et al. 2009) 
that w as found to be an AMXP with RXTE ( Markwardt et alj 
I2009ah . The INTEGRAL-deiived source position for the new 
transient was ~ 20' away from the known 435 Hz millisecond 
pulsar XTE J1751-305, the source position uncertainty of 2', 
as well as the subsequent detection of a distinctly different 
pulse frequency of 24 4.8 Hz confirmed the transient to be a 
new AMXP (Markwardt et all |2009a). The orbital period of 
IGR J17511-3 057 is ~ 3.47 hr, wifli an asin(0/c value of 
~ 275.2 It-ms dRiggio et al.ll201 ^. and references ther ein). A 
candidate near-infrared counterpart was found by Torr es et al] 
(l200 9a) on September 22, 2009 with a Ks-hand magnitude of 
18.0+0.1, which had faded lo K, > 1 8.8 (3cr upper limit) by 
October 7, 2009 dTorres et al.1 l2009bh . Radio upper limits of 
0.16-0.18 mJy between Se ptember 16 and 25, 2009 were set 
with the VLA by Miller- Jo nes et all (l2009h . The most accurate 
position of the source was provided by the near-infrared obser- 
yations at gponn = 17h51™08!64 and dmoo = -30°57'40'.'70 
dTorres et alJl2009al) . con sistent with the Icr error of 0'.'6 of the 
Chandm/UETG position (Nowa k et a0200 9^. 

T he first ty pe-I X-ray burst from the source was detected by 
Swift dBozzo et al.i 2009), an d the Swift m onitoring of this source 
in outburst was reported by iBozzo et"an 12010 ). who found no 
evidence for photospheric radius expansion during the three 
studied bursts. All follow-up observations with different space- 
crafts detected X-ray bursts from the source. Burst oscillations 
at ~ 245 Hz were first de tected in an X-ra y burst on September 
14, 2009 by RXTEfPCA dWatts et alJl2009i) : XMM-Newton also 



detected X-ray bursts wi th burst oscillation s phase locked to 
the persistent pulsations (Papitto et alJ l2009l) . A Chandra ob- 
servation on September 22, 2009 also detected an X-ray burst 
dNowak et alJ 2009). The net persistent spectrum, 0.5-10 keV 
band, of IGR J 1751 1-3057 was found for all the observations to 
be well-fit by an absorbed power-law model with photon index 
~ 2 and column density A^h - (0.6 - 1.3) x lO^^cm"^, depend- 
ing on the model u sed dBozzo et al J 12009: Nowak et al. 20091; 
iPapittoet alll2010ah . 

A lower hmit on the co mpanion mass of 0.1 3M0 was de- 
termined bv iMarkwardt et a l. (2009a) using the mass function 
of the system and assuming a NS mass of 1.4 M© and an in- 
clination angle of 90 °. This lower limit was later improved by 
iPapitto et al.l d2010al) . who considered that no occultation or dips 
were observed during the outburst of IGR J1751 1-3057. These 
authors also performed a fit to the source spectrum with a reflec- 
tion continuum model and determined an inclination angle for 
the system in the range 38 - 68°. The corresponding limits on 
the companion mass in this case are 0.15- 0.23 Mp (a NS of 1.4 
Mo was assumed). We refer the reader to 'Papi tto et al.l d2010ah 
for further detail regarding the nature of the companion star. An 
upper li mit on the source di stance of ~ 10 kpc w as estimated 
first by iPapitto et al.l (1201 Oal) and ~ 7 kpc by Altamirano et al.l 
d2010l) assuming that the most energetic type-I X-ray burst ob- 
served from IGR J1751 1-3057 reached the Eddington luminos- 
ity. 

In this paper we report on INTEGRAL and RXTE observa- 
tions of IGR J1751 1-3057, as well as a simultaneous Swift, 
INTEGRAL, and RXTE observation during the period MJD 
55087-551 17 (September 13 - October 13, 2009). We study the 
light curves, broad-band spectra, outburst spectral evolution, and 
timing properties of the source. The properties of the largest set 
of X-ray bursts from this source are also investigated. 



2. Observations and data 

2.1. INTEGRAL 

The present data were obtained during the INTEGRAL 
dWinkler et alJ l2003l) Target of Opportunity (ToO) observa- 
tion during satellite revolution 846, starting on September 16 
(MJD 55090.92557) and ending on September 19, 2009 (MJD 
55093.58993), with a total net exposure time of 206 ks. The 
data reduction was performed using the standard Offline Science 
Analysis (OSAQ) version 8.0 di stributed by the Integral Science 
Data Center ( Courvoisier et al ] l2b03i) . The algorithms u sed for 
the sp atial and spectral analysis are described in Goldwur m et al.l 
d2003h . The observation, aimed at IGR J1751 1-3057, consists 
of 70 stable pointings with a source position offset < 7^0 
from the center of the field of view. We analyzed data from 
the IBIS/ISGRI coded mask telescope dUbertini et all 120031; 
Lebrun et al. 2003) at energies betw een 18 and 300 keV and 
from the JEM-X monitor, module 1 dLund et al.ll2003l) between 
3 and 20 keV. Note, since this observation was in a rectangular 
pattern mode which consists of a square 5x5 pattern around the 
nominal target location, IGR J1751 1-3057 was only within the 
source position offset of < 3°5 from the center of the field of 
view for 28 stable pointings. Therefore, for JEM-X with a field 
of view of 7°5 (diameter) at half response the effective exposure 
time was only 87.4 ks. 

We first deconvolved and analyzed separately the 70 single 
pointings and then combined them into a total mosaic image in 



http://www.isdc.unige.ch/integral/analysis 
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the 20^0 keV and 40-100 keV energy bands, respectively. In 
the mosaic, IGR J17511-3057 is clearly detected at a signifi- 
cance level of ~ 65 cr for energies between 20^0 keV and ~ 29cr 
at higher energies (40-100 keV). The source position in the 20- 
40 keV band is q'j2ooo = 17f'51™08^.74 and 6moo = - 30°57'36V9 
(error of 0^55 at the 90 per cent confidence level, iGros et al] 
I200 3''). which is offset with respect to the near-infrared position 
by 0.'07 (Torres et al. 2009a). 

Figure [U shows part of the ISGRI field of view (significance 
map) around the position of IGR J 175 1 1-3057 in the 20-70 keV 
energy range. IGR J1751 1-3057 is clearly detected, together 
with two other nearby sources, at a high significance level (white 
circles). In the same figure we also indicated with a yellow cir- 
cle the position of the millisecond pulsar XTE J 175 1-305 that 
underwent a short period of enhanced X-ray activity during the 
last part of the outburst decay of IGR J1751 1-3057 (see Sect. 
IZHandFig.ia. 

To search for X-ray bursts, the ISGRI light curves are calcu- 
lated from events selected according to the detector illumination 
pattern for IGR J1751 1-3057. For ISGRI we used an illumina- 
tion factor threshold of 0.6 for the energy range 18^0 keV; for 
JEM-X we used the event list of the whole detector in the 3-20 
keV energy band. 

2.2. RXTE 

We used publicly availab le data from t he pro portional counter 
array PCA (2-60 keV; iJahoda et all 41996!)) and the High 
Energy X-ray T iming Experiment HEXTE (15-250 keV; 
[Rothschi idetaP (Il998l) ) on-board the RXTE satellite. The 
FWHM of these collimator instruments is ~ 1°, and no 
spatial information of the photons exists. IGR J175 11-3057 
was monitored from September 13 to October 8, 2009 (MJD 
55087.85619-55112.33693) for a total net exposure time of 
~ 455 ks (observation ID 94041). For the light curve analy- 
sis, we also added the twelve pointings from observation ID 
94042 (~ 45.2 ks net exposure time; target XTE J1751-305). 
In this dataset the nearby X-ray millisecond pulsar XTE J 1751- 
305 underwent a short outburst contemporaneous with the lat- 
est outburst phases of IGR J1751 1-3057 (Chenevez et al. 2009; 
iMarkwardt et al.l l2009b). Figure [1] shows that the angular sepa- 
ration between the two sources is smaller than the field of view 
of the non-imaging instruments on-board RXTE (~ 1°), and thus 
we were unable to disentangle the contributions of the two ob- 
jects to the total X-ray flux. 

We carried out a spectral analysis, using Standard-2 data 
(with 16 s time resolution) for the PCA and standard Cluster-0 
data for HEXTE. For HEXTE, we used the on-source data, using 
default screening criteria for Cluster 0. The data were extracted 
for the Good Time Intervals defined by standard criteria. The 
PCA response matrix was created by FTOOLS version 6.0.2 for 
129 energy channels to cover the energy range from 2-60 keV. 
For HEXTE we used the standard 64 energy channel response 
matrix for the 15-250 keV energy range. For the timing analy- 
sis we used PCU data, collected in the E_125us_64M_0_ls event 
mode, recording event arrival times with 122.07 fis time resolu- 
tion, and sorting events in 64 PHA channels. Default selection 
criteria were applied. 

2.3. Swift 



1 [depree] 

Fig. 1. INTEGRAL/lSGRl sky image of the field of view around 
IGR J1751 1-3057 (20-70 keV). The size of each pixel in the 
image corresponds to 3'. It is evident from the image that the 
non-imaging instruments on-board RXTE with a field of view of 
~ 1° (radius) would not be able to separate the emission from 
the two AMXPs XTE J1751-305 and IGR J1751 1-3057 when 
they are both active in X-rays (see Fig.|2]l. 



IGR J 1751 1-3057 was also observed with Swift dBurrows et al.l 
I2005h during its outburst in 2009. The results were first pub- 



Ushed in iBozzo et alj ( l2010h . We used only Swift/XRT data 
in window-timing mode (WT) that were simultaneously col- 
lected with the INTEGRAL observations (i.e. Swift observation 
ID 00031492005; see Table J_ in Bozzo et al. 2010, for details). 
We refer the reader to lBozzo et al.1 (l2010l) for the Swift XRT data 
reduction procedure. 

3. Outburst light curve 

Most AMXPs that underwent an outburst for a few weeks to 
months showed a common outburst profile, i.e. the light curve 
decays exponentially until it reaches a break, afte r which the flux 
drops Unearly to the quiescence level (see e.g., Gilfanov et al.l 
ll998HGierUnski & Poutanenll2005LlFaranga et al.ll2005bl) . These 
outburst profiles have been modeled fo r soft X-ray tr a nsient s 
based on the disk instability picture of King & Ritted (Il998l) . 
i.e. taking into account the disk irradiation by the central X-ray 
source during the outburst. iKing & RitteJ (Il998h showed that 
X-ray heating during the decay from outburst causes the light 
curves of transient LMXBs to exhibit either exponential or lin- 
ear declines depending on whether or not the luminosity is suf- 
ficient to keep the outer disk edge hot. This model has been ap- 

S >lied to a sarnple o f different LMXBs, inc luding three AMXPs 
Powell et al.]l2007h . iHartman et al.1 (l2010h investigated the out- 
burst light curve for one of the AMXPs, IGR J0029 1 +5934, us- 
ing different mo dels in which the d ecay tail is not necessarily 
linear. However, iPowell et al.l (l2007h showed that the timescale 
of the decay light curve and its luminosity at a characteristic time 
are linked to the outer radius of the accretion disk. A knee in 
the light curve of the decay from outburst of a transient LMXB 
is believed to be a consequence of mass transfer onto the outer 
edge of the disk, since this supply is effectively cut off from the 
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Fig. 2. RXTE/PCA (2-20 keV) outburst light curve of 
IGR J17511-3057. For plotting purpose we chose a bin time of 
6840 s. The count rate has been converted into flux using the 
spectral results reported in Sect. 14.21 The diamonds correspond 
to observations in which both XTE J 175 1-305 and IGR J 175 1 1- 
3057 were active and the instruments on-board RXTE were un- 
able to disentangle the contribution of the two sources. The ar- 
rows indicate the times of the detected X-ray bursts (see Sect. 
|6]and TableO. The vertical dashed lines indicate the interval of 
the INTEGRAL observations. The dot-dashed lines correspond to 
the best-fit exponential profile, F oc e '^22 4''^ ^jj^j ^ linear decay. 



compact object when the outer disk enters the cool low-viscosity 
state. The X-ray luminosity at which the knee occurs is that at 
which the outer disk edge is just kept hot by central illumination, 
allowing this radius to be calculated. In addition, the exponen- 
tial time-scale of the decay gives a second measure of the disk 
radius. 

The 2-20 keV band RXTE light curve of IGR J1751 1-3057 
has been extracted for all the pointings and is shown in Fig. |2l 
averaged over 6840 s. The diamonds correspond to the sum of 
the emission of the X-ray millisecond pulsar XTE J1751-305 
in a weak outburst and of IGR J1751 1-3057 returning back to 
quiescence. RXTE and INTEGRAL detected ten and three type-I 
X-ray bursts, respectively. These are indicated with arrows and 
are subtracted from the light curve and discussed in Section |6l 
We converted the count rates to flux using the spectral results 
from Section|421 The dot-dashed lines correspond to the best-fit 
exponential profile F oc g '/22 4 ^jj^j ^ linear decay. In order to 
fit the light curve of the outburst decay of IGR J 1751 1-3057, we 
excluded the RXTE data in which the other AMXP XTE J1751- 
305 was also active. We have verified that a linear fit would be 
probably favoured by the data (the reduced is about half of 
the coiTesponding value for the exponential fit). However, given 
the relatively poor observational coverage over the last part of 
the outburst decay, a firm distinction between a linear and an 

exponential decay cannot be made. 

From the equations (9), (23), and (21) in lPowell etal] (l2007h 
it is possible to estimate the outer disk radius in two difl'erent and 
independent ways. Equation (9) gives Lx - (L, - ) exp(-(f - 
tt)lTe) + Le, where L^, tt (the break time), L,, and t^. (exponen- 
tial decay time) are all free parameters. By using for these pa- 
rameters the values determined from the fit to the light curve of 
the outburst decay of IGR J1751 1-3057 (see Fig. the outer 
disk radius can be estimated from Rdiskije) - ije 3vkr)'^^ ~ 
4.8 X lO'" cm. Here we adopted for the viscosity near the 
outer disk edge, vtcR — 4 x 10'"^ cm^ s~', i n agreement with 
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lKing&RitteJ(ll998h and lPowell eTaP (l2007h . An independent 



Fig. 3. Unfolded broad-band spectrum of IGR J 175 11-3057 fit 
with the thermal Comptonization model compps. The data points 
are from XRT (0.8-7 keV), PCA (3-22 keV), and ISGRI (20-300 
ke V) spectra, respectively. The total model spectrum is shown by 
a solid line. The lower panel shows the residual between the data 
and the model. 



estimate of R jisk can also be obtained by using the equation 
(21) in Po well et all (l2007h . \.t RdiskiL,) = ('^L,y '^. By assum- 
ing O =s 1.3 cm^ s erg"' (see iPowell et al.ir2007h . we obtained 
RdiskiLi) ~ (4.8 - 5.4) X 10^" cm, for a distance to the source of 
6.3 and 7 kpc, respectively (see also Sect.|6]l. Note that the two 
independent estimates of the inner disk radius agree remarkably 
well in case a distance to the source of 6.3 kpc is considered. 
A similar ag reement was a lso fo und for a number of different 
AMXPs, see iPowell et al.1 (l2007l) . This radius also fulfills the 
condition Rd^ < Rdkk < b\, where R^m ~(2.8-2.3)xl0"' 
cm is the circularization radius and b\ ^ (6.7 - 6.4) x 10'" cm 
is the distance of the Lagrange point L \ from the center of the 
neutron star (see e.g.. iFrank et al.ll2002l) . To estimate these val- 
ues we used a companion star mass between 0.15 - 0.23 Mq 
(Pa pitto et al.,.2010a) . for a NS mass of 1 .4 Mq. 



4. Spectral analysis 

The spectral analys is was carried out using XSPEC version 
12.6 (lArnaudlll996h . For the contemporaneous data we com- 
bined the 3-22 keV RXTE/PCA data, and the 5-300 keV 
INTEGRAL/JEM-XflSGRl data taken on September 16-19, 
2009. In addition, for the low energy range 0.5-10 keV we also 
used the contemporaneo us Swift data (ID 00031492005 WT) 
from iBozzo etal.l(l20ltf ). For each instrument, a multiplication 
factor was included in the fit to take into account the uncertainty 
in the cross-calibration of the instruments. For all the fits the 
factor was fixed at 1 for the ISGRI data. To follow the outburst 
spectra outside the INTEGRAL observation interval we used the 
RXTE /PC A/HEXTE data. To take into account contamination 
by the Galactic ridge emission, we used as background for the 
PCA spectral fitting the RXTE/PCA data collected when both 
pulsars were at the lowest X-ray emission level (6.3 cts/s) (see 
also Papi tto et a02010atlAltamirano et al.ll2010l) . All uncertain- 
ties in the spectral parameters are given at a 90% confidence 
level for single parameter. 



4. 1 . Broad-band spectrum 

We studied in detail the broad-band X-ray spectrum of 
IGR J1751 1-3057 in the energy range 0.8-300 keV, us- 
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ing the joint INTEGRAL/JEM-X/lSGRl, RXTE/PCA, and 
Swift/XRT data. For the RXTE data we removed the time 
intervals corresponding to the bursts. We first fit the joint 
XRT/JEM-X/PCA/ISGRI (0.8-300 keV) spectrum using a sim- 
ple photoelectrically-absorbed power-law, pl, model which was 
found to be inadequate with a;t'^/d.o.f. = 417/141. A better 
fit was obtained by replacing the pl with the exponentially cut- 
off power-law model, cutoffpl, (x'/d.o.f. = 115/140). We es- 
timated a power-law photon index of F = 1.61 + 0.04 and a 
cutoff energy of = 58 ± 7 keV. The absorption was found to 
be A^H - 0.96!Q,)g x lO^^c m"^, close to the Galactic value re- 
ported in the radio maps of Dic key & Lockmaiil (1199 0'). Adding 
a multi-temperature disk blackbody model, diskbb, to the fit in 
order to include the contributio n of the soft therm al disk emis- 
sion to the total X-ray flux (iMitsuda et al.l [198 41) gave only a 
marginal improvement to the fit (x^/d.o.f. = 105/146). For this 
spectral component we obtained an inner disk temperature of 
kTin = 0.13 ± 0.02 and an inner radius of Ri„ ~ 601 Vcos ; km , 
in agreement with the results reported bv iBozzo et al.l (l2010h . 
The best fit value of Rin suggests a relatively large value of the 
inner disk radius that is hardly compatible with the idea that in 
AMXP sources the accretion disk extends down to a region very 
close to the NS surface (i.e., 10 km). A better description of the 
emission from IGR J 1751 1-3057 at the lower energies (<2 keV) 
was obtained by using the higher spectral resoluti on of the in- 
strumen ts on-ho3ird XMM-Newton, as discussed in iPapitto et alJ 
(I2010ah . 

In order to compare the IGR J175 11-3057 spectrum with 
previously observed broad-band spectra of the same source 
class (e.g., [O ierliriski et al. 2002; Gierlihski & Poutanen 2005; 
iFalanga et aP i2005a>ibi. i2007; : Jbragimov & Poutanen i2009i) . 
we replaced the cutoff power-law model with the ther- 
mal Comptonization model , compps, in the slab geometry 
dPoutanen & Svenssonlll996l) . The main model parameters are 
the Thomson optical depth Tt across the slab, the electron tem- 
perature kTfi, the temperature kTb\, of the soft seed blackbody 
photons assumed to be injected from the bottom of the slab, the 
emission area Abb, and the inclination angle between the slab 
normal and the line of sight. The best fit parameters of all the 
models used to fit the data are reported in Table [1] The marginal 
discrepancy in flux is due to the difference in the absorption col- 
umn density measured from the fits with the two spectral models. 
In the following we use A^h = 0.6 x lO^^cm"^, as obtained from 
the fit with the compps model. In Fig. [3] we show the unfolded 
spectrum and the best fit compps model. The residuals from the 
fit are also shown. In this fit, the normalization of the ISGRI data 
was fixed to 1 as a reference, while the normalizations of the 
XRT and PCA data were 1.02 + 0.08 and 1.05 ± 0.05, respec- 
tively. A reasonable value for the normalization of the RXTE 
PCA could be obtained only once the Galactic ridge emission 
was taken into account in the background correction of this data. 



4.2. Spectral evolution during outburst 

We analyzed all 58 RXTE /PC A/HEXTE (3-200 keV) spectra ob- 
served during the outburst. We excluded the last twelve RXTE 
pointings that were contaminated with the emission from a short 
weak outburst of XTE J1751-305 (see Fig. [T] |2l and Sect. 
O. In Fig. |4] we show the best-fit results using the thermal 
Comptonization compps model. The burst intervals have been re- 
moved. For the compps model the absorption column density and 
the inclination angle, 6, were fixed at the best-fit value found for 
the broad-band fit (see Table[TJ. The luminosity Lboi was calcu- 
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Fig. 4. The outburst evolution of the best-fit spectral parameters 
of the compps model using the RXTE /PC A/HEXTE data in the 
2.5-200 keV range. Each spectrum corresponds to one RXTE 
observation of IGR J1751 1-3057. The bolometric luminosity is 
given in the energy range 0.5-200 keV assuming a source dis- 
tance of 7 kpc. 

Table 1. Best-fit spectral parameters of the cutoff power-law and 
COMPPS models to the XRT/PCA/ISGRI data 





CUTOFFPL 


COMPPS 




96+" '" 


6+" ''' 

"•"-0.16 

24.8+^^ 


kT^ (keV) 




fcrtb (keV) 




59+"-°'* 


Tt 




9 9+0.(4 
^- -0.17 


Abb" (km2) 




260+**" 


cos 




53+" "'' 


r 


1 60+" 

-0.04 




E, (keV) 


58.3+^f 




^-/dof 


112/140 


106/138 


Fboi*(10-'' ergcm-2 s"') 


1.41+0.10 


1.28+0.11 



Notes. Assuming a source distance of 7 kpc. Unabsorbed flux in 
the 0.8-300 keV energy range. 



lated for a distance of 7 kpc (see Sect.|6]l from the flux measured 
for the best-fit model in the energy range 0.5-200 keV. These 
results show that the decay of the outburst is marked by a nearly 
constant plasma temperature, feT^e, soft seed photons emission, 
kT\,b, and optical depth, t. No statistically significant variations 
are measured. 

A similar behaviour, althrough less evident due to less ob- 
servational coverage during the outburst, was also observed i n 
the case of the AMXP IGR J0029 1+5934 (Falanga e t al]l2005ah . 
Note that IGR J1751 1-3057 and IGR J0029 1 +5934 show simi- 
lar co l or variations during the outburst decay dAltamirano et al.l 
120101 : iLinares et al.l l2007h . The hardness-intensity diagram of 
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IGR J 175 11-3057 is shown in Fig. 2 of lAltamirano 6110(1201 Oh . 
and shows a stable (hard) color up to ~ 18 days after the on- 
set of the outburst, then the color softens within the hard state. 
Our Fig. |4] also shows that the product t x feT^e, as measured 
during the outburst of IGR J1751 1-3057 is stable during the 
outburst, in agreement with observations of other AMXPs (e.g.. 



Gilfanov et al. ' l998';'Gierliriski & Poutanen"2005';'Falanga et alj 
2005a; Poutanen 2006; Ibragimov & Poutanen 2009). Such be- 
haviour is expected if the energy dissipation takes place in 
an accretion shock. The shock geometry can be approximated 
by a slab and the cooling of the hot electron gas is deter- 
mined by the reprocessi ng of the hard X-ra y radiatio n at the 
neutron star surface (see Haar dt & Maraschi 1993; Stern et alJ 
ll995l;lPoutanen & Svensson 1996; Malzac et alj2001h . The tem- 
perature depends on the optical depth, but tj x kT^ is ap- 
proximately constant. In IGR J1751 1-3057 we observe kT^ x 
tt ~ 50 keV (see Fig. |4]i whic h is consistent with the val- 
ues deter mined for other AMXPs (Poutanen & Gierlinskill2003l: 
Gierliris ki & PoutanenI llool iFalanga et al.l l2005allbl ; iPoutanenI 
2006) as well as with the theoretical models. 



The spectra of AMXPs are very simil ar to th ose of the atoll 
sources at low luminosities (iBarret et all |200() I). where the X- 
rays are probably produced in the boundar y/spreading layer near 
the NS equator ( Kluzniak & Wilson 1991; Inogamov & SunvaevI 
Il999t ISuleimanov & Poutanenl i2006). Spectral similarities can 
be explained if in both types of sources, the energy dissipation 
happens in the optically thin medium (i.e. accretion shock and 
boundary/spreading layer) and the spectral properties are deter- 
mined solely by energy balance and feedback from the NS sur- 
face which provides cooling in the form of soft photons. 



5. Timing characteristics 

We also studied the pulsed emission of IGR J1751 1-3057 in 
the 3-300 keV energy range using RXTE/PCA, HEXTE, and 
INTEGRAL/lSGRl data. As the flux decays over the course of 
the INTEGRAL observation we selected only those RXTE ob- 
servations that overlap with the INTEGRAL observation lasting 
from MJD 55090.926 to 55093.520 (UTC) (see Fig.|2|. In addi- 
tion to the application of standard selectioncriteria in the screen- 
ing process the RXTE/PCA data were further screened for bursts 
and detector break-downs. The screening process yielded the 
following exposure times for PCU detectors 0^, respectively: 
16.096 ks, 12.960 ks, 47.344 ks, 10.848 ks, and 13.216 ks. For 
a discussion of ti ming behaviour over t h e whole o utburst we 
refer the reader to Pa pitto et aP (1201 Oah : iRiggio et alJ (12010l) : 
[ibragimov et a l. (2pifl). 

For HEXTE we only selected the on-source data streams 
from both detector clusters. The dead time corrected exposure 
times in this case are 30.2 and 16.0 ks for HEXT^I Cluster and 
1, respectively, for the period overlapping with the INTEGRAL 
observation. For the full 94041 observation the effective expo- 
sure times were 232.4 and 123.2 ks for Cluster and 1, respec- 
tively. 

We screened the ISGRI data for effects due to Earth radiation 
belt passages or solar flare activity, none of the pointings showed 
such activity. Time periods in which burst events occured from 
any source in the ISGRI field of view were excluded from further 
analysis. Furthermore, we selected only time stamps of events 
with rise times between channels 7 and 90 (Lebrun et al. 2003), 
detected in non-noisy detector pixels which have an illumination 
factor of more than 25%. 
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HEXTE Cluster operates in staring mode since July 13, 2006 



Fig.5. Pulse-profile coflage of IGR J17511-3057 using data 
from RXTEfPCA, HEXTE, and INTEGRALflSGRl. Panels a-d 
(PCA) correspond to the energy intervals lJ-4.2, 4.2-7.9, 7.9- 
16.2 and 16.2-34.9 keV. Panels e-g show the HEXTE pulse pro- 
files in the bands 15.6-31.0, 31.0-60.1 and 60. 1-123.9 keV (ac- 
cumulated during the entire 94041 observation period) in panel h 
the ISGRI profile 15-90 keV. The error bars represent Icr statis- 
tical errors. All profiles reach their maximum near phase ~ 0.1. 
Note, the y-label are in units of counts per bin. 



The selected time stamps of all the instrument were then con- 
verted to arrival times at the Solar System barycentre taking into 
account the orbital motion of the spacecraft and correcting for 
acceleration effects along the binary orbit. In this process we 
used the positi on of the opt i cal cou nterpart to IGR J1751 1-3057 
as reported bv lTorres et all ( l2009allbl) . 



5.1. Pulse profiles and time lags 

Fhilse phase foldin g of the barycentred arrival times using the 
ephemeris given in lPapitto et al.l (l2009l) yielded the pulse-phase 
distributions shown in Fig. |5] In case of the PCA data (col- 
lected in mode E_125us_64M_0_ls) we added a time shift of 
+0.5 X 122.07yus to the barycentred time stamps because the 
times refered to the start of the time bin instead of the mid 
of the bin. Panels a-d of Fig. |5] show the PCA pulse profiles 
for the energy intervals, 1.7-4.2, 4.2-7.9, 7.9-16.2 and 16.2- 
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34.9 keV, respectively. In panels e-g the HEXTE profiles (ac- 
cumulated during the full 94041 observation) are shown for 
the bands 15.6-31.0, 31.0-60.1 and 60.1-123.9 keV, and panel 
h shows the ISGRI profile for energies between 15 and 90 
keV. Pulsed emission has been detected up to ~ 120 keV us- 
ing HEXTE data collected during the full duration of observa- 
tion 94041. The mutual alignment of the RXTE/PCA, HEXTE, 
and INTEGRALflSGRl profiles within the equivalent energy 
bands is better than 50jus. The profiles are rather sinusoidal with 
the amplitude s of harmonics not exceeding 15% of the fun- 
damental (see iPapitto et al.ll2010ah . A small asymmetry could 
be caused by a strong deviation of emis sion pattern from the 
blackbody ( Potit anen & Giei'linski 2003: iViironen & PoutanenI 
|2004|) and/ or the appearance at some phases of the sec- 
ondary spot (Papitto et al."2010a') as seen in S AX J1808.4-3658 
Jlbragimov & Poutanen 2009; Poutan en et"an i2b09). 

We studied the global arrival times of the pulses as a function 
of energy, and used the PCA profile for channel 1 (~ 1 .69 - 2.5 1 
keV) as a reference template. Cross-correlation of the pulse pro- 
files obtained for other energy bands with this template yielded 
the time lags shown in Fig.|6l Beyond ~ 4 keV a declining trend 
sets in, meaning that the hard X-ray photons arrive earlier than 
the soft ones. 

The behaviour of the IGR J175 11-3057 time lags is simi- 
lar to that observed in XTE J175 1-305. In these cases the low- 
energy pulses lag behind the high-energy pulses (soft phase/time 
lags) mono tonic all y increasing with energy and saturating at 
about 10-20 keV dOierliriski & PoutanenI l2005h . The satura- 
tion ene rgy is, howeve r, much smaller in IGR J0( )291-h 5934 
(see e.g. Falanga et al. 2005b; Falanga & Titarchuk 2007) and 
SAX J1808.4^3658 (Cui et al. 1998- Pout anen & Gierlihskl 
120031; lPoutanenll2006l: llbragimo v & PoutanenH2009l) . where the 
time lag as a function of energy breaks at about 7 keV. For 
IGR J0029 1-1-5934 the soft lag even decreases beyond. 

The lags are most probably due to different emis- 
sion patterns of the blackbody and Comptonization com- 
ponen ts (Poutanen & Gierlihski 200jl iGierlihski & PoutanenI 
120051; llbragimov & Poutaneni2009i) combined with the action of 
the Doppler effect. This is supported by the energy dependence = 
of the lags, which grows until the contribution of the blackbody 
becomes negligible. At higher energies, the evolution of the soft 
lag is still further possible, because higher-energy photons suffer 
more scatterings resulting in variations of the emission pattern 
with e nergy (Poutanen & Gierlihski 2003; Viironen & Poutane^ 
l200l . The break in the lag spectrum at a higher energy 
might be a result of a higher seed blackbody temperature for 
Comptonization, which in turn implies a smaller spot area, 
and/or a smaller contribution of the blackbody to the total 
spectr um. Note that in XTE J1751-305 (Gierlihski & Poutan^l 
I2OO5I) and IGR J1751 1-3057 the blackbody is hardly visible, 
while it is very appa rent in the spectra of IGR J0029 1-1-5934 
(Falang a^l.ll2005bh. and SAX J1808.4-3 658 ( GierHhskTetaD 
[20021; iPoutanen & Gierhhskil 120031; iPoutan en 2008). A smaller 
blackbody area indicates either a larger magnetic field in the 
first two sources compared to IGR J0029 1-1-5934, or a differ- 
ent geometry (e.g. inclination), or a different optical depth of the 
Comptonizing plasma in the accretion shock. 

5.2. Spectrum of the pulsed emission and pulsed fraction 

The pulsed spectrum and pulsed fraction (defined as pulsed 
flux/total flux) as a function of energy provide important diag- 
nostic parameters for cons training the parameter spac e in the- 
oretical modeling (see e.g. lViironen & Poutanenll2004l) . We de- 
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Fig. 6. Time lag as a function of energy in the ~2-100 
keV energy range combining RXTE/PCA (1.7-35 keV; open 
circles), RXTE/HEXTE (15-90 keV; filled diamonds), and 
INTEGRAL/lSGRl (15-90 keV; fifled squares) measurements. 
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Fig. 7. Unabsorbed total (-0.8-300 keV) and pulsed (-3-120 
keV) unfolded spectra of IGR J1751 1-3057 combining mea- 
surements from RXTE/PCA, HEXTE, and INTEGRAL/lSGRl 
for the pulsed part and Swift/XRT, INTEGRAL/JEM-X and 
ISGRI for the total part. 



rived the pulsed fluxes (~3-120 keV) from observations with 
the non-imaging RXTE/ PCA and HEXTE instruments overlap- 
ping in time with the INTEGRAL observation, and the (imaging) 
ISGRI instrument. The total fluxes (-0.8-300 keV) have been 
determined from instruments with imaging capabilities, namely 
SwiftfXRT, and JEM-X and ISGRI aboard INTEGRAL (see also 
Sect. El]). 
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Firstly, we derived the pulsed excess counts (=counts above 
DC level) in a given energy band by fitting a truncated Fourier 
series, using only the fundamental and two harmonics, to the 
measured pulse phase distribution. For the PCA data (~3-35 
keV band) the pulsed excess counts were converted to pho- 
ton flux values (ph cm^^s"' keV"') in a forward spectral fold- 
ing procedure assuming an underlying power-law model taking 
into account the different exposure times (see Sect. |5]l and en- 
ergy responses of the five active PCU's. We kept the absorbing 
hydrogen column density fixed to 6.0 x 10^' cm"^ (see Table 
[TJ, although its precise value has very little impact on the fit 
results because we consider only measurements with energies 
above 3 keV. We obtained an unabsorbed 2-10 keV (pulsed) 
flux of 8.09(6) X 10 " erg cm^^s ' and a power-law index of 
1.824 + 0.004 (reduced of the fit was 14.34/(21-2)), indicat- 
ing that the pulsed spectrum is softer than the total spectrum (see 
Table[B. 

For the HEXTE data (~ 15-120 keV) we employed an equiv- 
alent method: the pulsed excess counts in a certain energy band 
(for those RXTE observations overlapping with the INTEGRAL 
observation) were divided by its effective sensitive area assum- 
ing a power-law model with index 1.824, taking into account 
the different energy responses and dead-time corrected exposure 
times of the two detector clusters. Finally, the pulsed ISGRJ ex- 
cess counts (-15-90 keV) of IGR J1751 1-3057 have been con- 
ve rted to flux values ad opting the method outlined in Section 3.4 
of lKuiperetal] (120061) . 

The (unabsorbed) PCA, HEXTE and ISGRI pulsed flux mea- 
surements are shown in Fig.Qalong with the (unabsorbed) total 
flux measurements. If we compare the energy spectrum of the 
puls ed emission of IGR J 17511-3057 and IGR J00291+5934 
(see iFalanga et al.l 12005 al) . then we see a significant differ- 
ence: in the EFe spectral representation the pulsed flux of 
IGR J175 11-3057 shows a cut-off^ near 20 keV, while that of 
IGR J0029 1-1-5934 increases up to the end of the sensitivity 
window. From the pulsed and total flux measurements shown 
in Fig. |2]we can directly derive the pulsed fraction as a func- 
tion of energy which is shown in Fig. |8] The pulsed fraction 
decreases from ~22% at 3 keV to a constant pulsed fraction 
of ~17-18% between 7-30 keV, and then decreases again un- 
til ~13% at 60 keV. The decreasi ng part is very similar to wha t 
is observed in XTE J1751-305 (iGierUnski & Poutanenll2005l) . 
while in IGR J0029 1-1-5934 the pulsed fraction has a minimum 
of about ~6% at 7 keV and grows again to ~ 12-20% at 100 keV 
(IFalanga et al.ll2005ah . 

6. Properties of the X-ray bursts 

In Table|2]we report the key measurable parameters for the bursts 
observed from IGR J175 1 1-3057. Thermonuclear (type-I) X-ray 
bursts are produced by unstable burning of accreted matter on 
the NS surface. The spectrum from a few keV to higher ener- 
gies can usually be well described as a black-body with temper- 
ature kTbb ~l-3 keV. The energy-dependent decay time of these 
bursts is attributed to the cooUng of the NS photos phere result- 
ing in a gradual softening of the burst spectrum (see lLewin et al.l 
[l993; Strohmaver & Bildsten 2006, for a review). 

We defined the burst start time as the time at which the X-ray 
intensity of the source first exceeded 10% of the burst peak flux 
(above the persistent intensity level). The time-resolved spec- 
tral analysis of the 13 bursts (10 RXTE and 3 INTEGRAL) was 
can-ied out by using RXTE/PCA and INTEGRAL/]EM-X data in 
the 2.5-20 keV and 3-20 keV bands, respectively. From these 
analyses we determined the bursts' peak fluxes, temperatures. 
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Fig. 8. The pulsed fraction (=pulsed flux/total flux) of 
IGR J1751 1-3057 based on the pulsed/total flux measurements 
from Swift/XRT, RXTE /PC A/HEXTE, and INTEGRAL/iEM- 
X/ISGRI. 



and radii (see Table |2]l. We fitted each burst spectra by a sim- 
ple photoelectrically-absorbed black-body, bb, model. The neu- 
tral absorption column density A^h was left free to vary in the 
range (2.6-0.001)10^^ cm"- in all fits. However, we checked that 
fixing A^H - 0.6 X 10^^ cm"^ would not significantly affect the 
results. We extrapolated the unabsorbed fluxes to the 0.2-50 keV 
band by generating dummy responses (xspec version 12.6). This 
is justifiable for the data since the black-body temperature is well 
inside the spectral bandpass. The infeiTed bb peak temperature, 
^T^bb.peak, apparent bb radius at 7 kpc (see below), /?bb,peak, and 
unabsorbed bolometric peak flux are also reported in Table|2| In 
Fig. |9| we show representative time-resolved spectroscopic pa- 
rameters. 

We calculated the burst fluence, /b, by integrating the flux 
over the burst duration. The effective burst duration is t\, = 
/b/Fpeak- All bursts are short, with t^, in the range 7-10 s. The 
fluence and peak flux for the bursts increased steadily with time, 
from ~ 3.0 x 10"^ to 4.2 x 10"^ erg cm"^ s"' for the fluence, 
and from ~ 3.3 x 10"^ to 5.7 X 10"^ erg cm"^ s ' for the peak 
flux. As a consequence, Tt, decreased steadily, and the rise time 
also decreased, from ~ 2 to ~ Is. The brightest burst reached 
(5.7 + 0.12)xl0"** erg cm"^ s"', and the minimum observed sep- 
aration of the bursts was 7.08 hr (se e Table O. These result s 
are in agreement with those reported in lAltamirano et al"] ( l2010l) . 
Those authors also reported the properties of burst oscillations 
detected at the pulse frequency in all the bursts. 

When a burst undergoes a photospheric -radius expansion 
(PRE), the source distance can be determined based on the as- 
sumption that the bolometric peak luminosity is saturated at the 
Eddington Hmit, LEdd (e.g., iLewin et aDll993t iKuulkers et al.l 
2003). None of the observed bursts exhibited PRE; thus, in 
these cases the peak luminosity is expected to have been sub- 
Eddington, leading to an upper limit on the distance. Assuming 
a bolometric peak luminosity equal to the Eddington value for a 
He X-ray burst (Lndd ~ 3.8 x lO-''^ ergs"', as empirically derived 
by Kuulkers et al. 2003), we obtain, using the brightest burst 
(see Table |2]l, the source distance upper- limit of li < 7.5 kpc. 
For comparison, the theoretical value of this upper limit distance 
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Table 2. Burst characteristics observed with RXTE/PCA and INTEGRAL/JEM-X 
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Notes. Pre-burst unabsorbed flux in 0.8-200 keV energy range. Burst peak flux in 0.1-40 keV energy band. '''' Burst fluence in 0.1-40 
keV energy band. Effective duration t\, = /b/-Fpcak- '"^^ Burst peak temperature. Burst peak blackbody radius for distance of 7 kpc. Bursts 
detected with 1NTEGRAL/]EM-X. 
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Table 3. Bursts recurrence times and effective exposure times. 
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Fig. 9. Representative time-resolved spectroscopic results from 
the brightest burst seen by RXTE from IGR J1751 1-3057. 
Shown is the inferred bolometric luminosity, from integrating 
the best-fit blackbody model (top panel); the blackbody temper- 
ature [middle panel); and the effective blackbody radius (bottom 
panel), assuming a distance of 7 kpc. 



(e.g. JLewin et al ■ll993h foundbv assuming a He atmosphere and 
canonical NS parameters ( 1 .4 solar mass and radius of 1 km), is 
~ 6.3 kpc. Alternatively, assuming the peak luminosity to be the 
Eddington luminosity for solar composition (X = 0.7) implies 
a limit of ^ 4.8 kpc. In the following, we consider d ^1 kpc 
to be a fiducial distance. At this distance, all bursts occuiTed at 
persistent luminosities between (5.2 - 9.6) x 10^^ erg s"' (see 



Table IS, or ^ (1.4 - 2.5)%LEdd using Lndd ~ 3.8 x 10^** erg 
s The local accretion rate per unit area for the pre-burst emis- 
sion, Lpers, is then given by m - Lpers(l +z)(4;r/?^(GM//?))"', i.e. 



# 


Observatory 


^stait 


Af,ec (hr) 


T'exp (hr) 


A' 


1 


Swift 


Sep 14 00:50:27 








2 


RXTE 


Sep 14 07:54:43 


7.084 


1.515 




3 


RXTE/Swift 


Sep 15 17:17:23 


33.395 


10.454 


3 


4 


RXTE 


Sep 17 06:33:27 


37.253 


10.88 


3 


5 


RXTE 


Sep 17 14:48:42 


8.254 


6.549 




6 


INTEGRAL 


Sep 18 08:09:29 


17.331 


7.859 


1 


7 


INTEGRAL 


Sep 19 02:28:20 


18.314 


7.303 


1 


8 


INTEGRAL 


Sep 19 11:13:03 


8.745 


8.226 




9 


RXTE 


Sep 20 14:50:31 


27.639 


5.056 


2 


10 


XMM-Newton 


Sep 21 01:12:21 


10.349 


2.40 




11 


XMM-Newton 


Sep 21 12:34:01 


11.361 


11.361 




12 


Chandra 


Sep 22 12:54:59 


24.349 


12.555 


1 


13 


RXTE 


Sep 23 14:27:06 


25.553 


8.207 


1 


14 


RXTE 


Sep 24 18.39:13 


28.214 


10.976 


1 


15 


RXTE 


Sep 25 07:31:36 


12.862 


6.133 




16 


RXTE 


Sep 26 15:11:22 


31.659 


13.126 


1 


17 


RXTE 


Sep 27 06:57:21 


15.770 


4.985 




18 


Swift 


Sep 30 18:31:57 


83.576 


24.240 


3 



Notes. The measured recurrence time between observed bursts is Af,.ec. 
The effective exposure time, Texp, takes into account all the instruments 
observations. In Fig.[lO]the triangles corresponds to ht^^^liN + 1), with 
A' being the number of the expected missed bursts. 



m ~ (2.13 - 5.3) X lO-'g cm"^ s"'. We use here the gravitational 
redshift I + z - 1.31 for canonical NS mass, M - I.4M0, and 
radius, R = IQ km. 

The observed energies of the bursts allow us to estimate 
the ignition depths. The measured fluences of the bursts are 
/b - (3.2 - 4.2) X 10"^ erg cm"^, corresponding to a net burst en- 
ergy release Eburst = 4;rt/Vb = (1.9 - 2.5) x 10^'^ (d/l kpc)^ erg. 
The ignition depth is given by jign = -EburstCl + z){4-7TR^Qnuc)''^ , 
where the nuclear energy generated (assuming a mean hydrogen 
mass fraction at ign ition {X}) is Qnuc ~ 1 .6 H- 4{X} MeV/nucleon 
(iGallowav et al.ll20 04. and references therein), including losses 
due to neutrino emission following Fuiimoto et al. ( 1987). Most 
burst sources accrete a mix of hydrogen and helium, although 
the H-fraction at ignition may also be reduced by steady burn- 
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ing between the bursts, i.e. X < Xq. Steady H-burning proceeds 
via the hot-CNO cycle, which will exhaust the available hydro- 
gen in approximately 11 hr (Z/0.02)"'(Xo/0.7), where Z is the 
mass fraction of CNO nuclei, and Xn the mass fraction of hydro- 
gen in the accreted fuel (e.g., iGallowav et al.|[2004.) . As for all 
bursts observed from IGR J1751 1-3057, the bursts early in the 
outburst, when the recurrence time was as short as 7 hr, exhibit 
short profiles, characteristic of low H-fraction, and with corre- 
spondingly high a-values. In order to exhaust the accreted H in 
just 7 hr via steady burning, the accreted H-fraction must be sub- 
stantially below solar; we infer Xq < 0.44(Z/0.02). Naturally, it 
is also possible that the H-fraction is closer to the solar value, but 
that the CNO metallicity is elevated. For the inferred pure helium 
composition at ignition (i.e. (X) - 0) the column depth varies lit- 
tle from burst to burst, in the range yign = (1.2-1.6)xl0'^ g cm"^. 

Once steady H-burning exhausts the accreted hydrogen at 
the base of the fuel layer, this process no longer dominates 
heating in the layer, and ignitio n will occur in a pure He-layer 
(case 2 of lFuiimoto et al.iri98lh . As the accretion rate decreases 
through the decay of the outburst, we expect the burst recuiTence 
times to become longer, as is observed. We found that the recur- 
rence time increases roughly as (/^pei-s.boi) ' ' (Fig. [TOb . where 
(Fpers.boi) is the averaged persistent flux between the bursts. 
However, the expected increase in the burst recurrence time in 
the pure He-ignition regime as a function of decreasing m is 
m uch steeper. For exa mple, the curve for = 0.3 in Fig. 1 
of lGallowav & Gummin g (2006) falls off roughly as m"^^, sub- 
stantially steeper than measured in IGR J17511-3057. What is 
also puzzling is that the measured a-value (the ratio of persis- 
tent to burst luminosity) decreases with decreasing (Fpers.boi)- As 
the burst recurrence time drops, the fuel layer should become 
increasingly He-rich, so that a slight increase in a would be ex- 
pected. 

This analysis relies on unambiguous measurement of the 
burst recurrence time, for w hich we used all the bursts reported 
in the literature, wit h Swift (iBozzo et alJl2009D. XMM-N ewton 
dPapitto etalj|2010al) . and Chandra dNowaketal] 120091) in ad- 
dition the INTEGRAL and RXTE. In some cases, the effective 
exposure time on source from burst to burst was shorter then the 
expected recurrence time between bursts, so that we most likely 
miss one or more bursts between the observed one (see Table 
[3]). Therefore, we divided some measured times between bursts 
by an integer number, Afrec/(A^ +1) (see Table O, where is 
the number of missed bursts. We verified that at the expected 
burst times we always had an observational data gap. With the 
inclusion of observations for which there were no intervening 
data gaps (e.g. with XMM-Newton), we are confident that we 
have precisely inferred the recurrence time between each pair of 
bursts listed in Table [3] 

The discrepancy between the expected and observed varia- 
tion in recurrence time is substantial, and in contrast to the (oth- 
erwise rather similar) bursts observed at comparable acc retion 
rates in SAX J1808.4-3658 (iGaUowav & Cummin d 120061) . As 
is often argued, it is possible that the variation in m is not so 
large as is suggested by the decrease in Fpers.boi- If, as the X-ray 
intensity drops, we are seeing a decreasing fraction of the accre- 
tion energy in the X-ray band, the actual m could remain higher, 
which would be consistent with both the moderate recurrence 
times and the decreasing a values. However, this would require a 
large change in the efficiency of the m to Fpers.boi conversion. One 
way this could occur is that if the spectral energy distribution 
changes such that the accretion flux is increasingly emitted out- 
side the band we are sensitive to. However, the combination of 
instruments with both low- and high-energy sensitivity (Swift & 
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Fig. 10. The triangles in the figure represents the observed burst 
recurrence times (see TableO shown as a function of local mass 
accretion rate (the corresponding flux is reported in the upper x- 
axis). We also show the best fit power-law model. The recurrence 
time increased with time roughly as (i^pers.boi) ' 



INTEGRAL) makes this unlikelv. IThompson et al] (l2008h found 
from a detailed study of GS 1826-24 that the characteristic un- 
certainty that might arise from this mechanism is 40%. This is 
not sufficient to explain the discrepancy in IGR J1751 1-3057. 

7. Summary 

We analyzed the simultaneous INTEGRAL, RXTE, and Swift ob- 
servations in order to study the broad-band spectrum and timing 
behaviour of IGR J1751 1-3057. Using all RXTE data we also 
studied the outburst profile. The broad-band average spectrum is 
well described by thermal Comptonization with an electron tem- 
perature of ~ 25 keV and Thomson optical depth rx ~ 2 in a slab 
geometry. The object shows remarkable spectral stability during 
the outburst marked by constant plasma and seed photon temper- 
ature at a constant scattering optical depth. We fitted the outburst 
profile with the exponential model and using the disk instability 
model we inferred the outer disk radius to be (4.8 - 5.4) x 10'° 
cm. 

We showed that the coherent pulsation can be tracked with 
HEXTE and ISGRI instruments up to ~ 100 keV. The pulsed 
fraction is shown to decrease from ~ 22% at 3 keV to a con- 
stant of ~17-18% between 7-30 keV, and then to possibly de- 
creasing further down to ~13% at 60 keV. This is similar to 
that observed in XTE J1751-305, but markedly different from 
IGR J00291-H5934. The nearly sinusoidal pulses show soft lags 
monotonically increasing with energy to about 0.2 ms at 10-20 
keV, with some indications of a further decrease. 

Using all observations by INTEGRAL, RXTE, Swift, 
Chandra, and XMM-Newton we have collected the largest set 
of X-ray bursts observed from IGR J1751 1-3057, which al- 
lowed us to determine the recuiTence time (accounting for the 
missed bursts) as a function of the accretion rate and the ignition 
depth. The short burst profiles indicate hydrogen-poor material 
at ignition, which suggests either that the accreted material is 
hydrogen-deficient, or that the CNO metallicity is up to a factor 
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of 2 times solar. However, the variation of burst recurrence time 
as a function of ih (inferred from the X-ray flux) is much smaller 
than predicted by helium-ignition models. 
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